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Abstract 

In this paper we consider the sequence of Kakutani's a-refinements 
corresponding to the inverse of golden ratio (which we call Kakutani- 
Fibonacci sequence of partitions) and associate to it an ergodic interval 
exchange (which we call Kakutani-Fibonacci transformation) using the 
"cutting-stacking" technique. We prove that the orbit of the origin 
under this map coincides with a low discrepancy sequence (which we 
call Kakutani-Fibonacci sequence of points), which has been also con- 
sidered by other authors. 
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1 Introduction and preliminaries 

In this paper we study with a dynamical system approach a sequence of points 
which arises from a beautiful geometric construction due to Kakutani [T2] : 
the a-splitting procedure. He proved the uniform distribution of successive 
a-refinements using ergodic theory methods. It has been proved in [2] that 
when a is the inverse of the golden ratio, the sequence of a-refinements 
(called Kakutani-Fibonacci sequence of partitions for reasons which will be 
clear later) has low discrepancy. 

In [2] and [3] the first author proposed two algorithms which associate 
to this sequence of a-refinements a sequence of points (called Kakutani- 
Fibonacci sequence of points) which has low discrepancy, too. 



In this section we present some basic definitions on uniform distribution, 
discrepancy, ergodic theory and the well-known Birkhoff 's Theorem. We also 
recall the definition of the Kakutani- Fibonacci sequence, with the addition of 
some properties, the most important of which being Theorem [6j 

In Section 2 we introduce an interval exchange T associated to the cutting- 
stacking procedure for this sequence and we prove that this transformation 
is ergodic (Theorem H%|) . By means of this result and Bikhoff's Theorem, we 
conclude that the orbit of almost every x G [0, 1[ is a uniformly distributed 
sequence of points (Theorem [T5l). With the last results (Theorem [16]) we 
obtain the Kakutani-Fibonacci sequence as the orbit of the origin under 
T. In a way, our results parallel Lambert's observation concerning the von 
Neumann-Kakutani mapping and the van der Corput sequence [To] . 

Both sections are endowed by figures which should be of some help for 
the reader. 

Definition 1. Let a G ]0, 1[ be a real number and n be any partition of [0, 1[ . 
Kakutani's a-refinement of 7r, denoted by air, is obtained by splitting all the 
intervals of 7r having maximal length into subintervals of lengths proportional 
to a and 1 — a, respectively. 

The sequence of partitions {a n u} n >i obtained by successive a- refinements of 
the trivial partition u — {[0, 1[} is called the Kakutani a-sequence. 

Definition 2. Let {7i n } n >i be a sequence of partitions of [0, 1[, represented 
by 7r n = {[yp 1 ,Vj^i[ : 1 < j < t n ). The discrepancy of {7r n } n >i is defined by 



We say that {n n } n >i is uniformly distributed (u.d.) if D{ii n ) — > when 
n — )• oo. 

If there exists a constant C > such that t n D(n n ) < C for any n, we say 
that {7r ra } n >i has low discrepancy. 

Definition 3. The discrepancy of the first N points of a sequence X = 
{x n }n>i in [0, 1[ is defined by 



D(7r n ) = sup - 

0<a<6<l t 



-^2xia,b[(y\ n) ) -ib-a) 



D N (X) = sup 

0<a<6<l 
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We say that X is uniformly distributed (u.d.) if D N (X) — > when n — > oo. 
If there exists a constant C > such that ND^(X) < C \ogN for any N, 
we say that X has low discrepancy. 

The following theorem shows the important role of u.d. sequences of 
points in Quasi-Monte Carlo methods. 

Theorem 4. A sequence of points {x n } n >i in [0, 1[ is uniformly distributed 
(u.d.) if and only if for any continuous function / on [0, 1] it is 



Kakutani proved that the sequence of partitions {a n u} n >i is uniformly 
distributed p2]. 

This result got a considerable attention in the subsequent decades and in 
the recent years many extensions of Kakutani's idea have been considered. 

The concept of u.d. sequences of partitions has been extended to compact 
metric probability spaces [7] and to fractals [H]. Recently, the first and third 
authors showed how to rearrange a dense sequence of partitions in order to 
obtain a u.d. one [6]. 

The third author [T7j extended Kakutani's a-refinement by introducing 
the concept of p-refinement of a partition n of [0, 1[. It is obtained by splitting 
the longest intervals of ir into a finite number of parts homothetically to a 
fixed finite partition p of [0, 1[. It has been proved that the sequence of 
subsequent /^-refinements {p n o>}„>i is u.d. for every p. 

Kakutani's splitting procedure has also been extended to the multidimen- 
sional case in [5]. 

The general case in which the successive p-refinements are applied to an 
arbitrary non trivial partition ir has been studied in [TJ. The authors give 
necessary and sufficient conditions on n and p which assure that {p n ^} n >i is 



Drmota and Infusino [8] derived bounds for the discrepancy of the se- 
quences of p-refinements {p n u} n >i considering a new approach based on the 
Khodak algorithm for parsing trees. 

A specific countable class of sequences of p-refinement s, called LS'-sequen- 
ces of partitions, has been studied in [2]. In that paper it was also presented 
an explicit algorithm which orders the points determining the LS'-sequences 




u.d.. 
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of partitions to obtain a u.d. sequence of points, called LS'-sequence of points, 
providing estimates for the discrepancy. 

The paper [3] contains a new explicit algorithm "a la van der Corput" 
to generate these sequences of points using the representation of natural 
numbers in base L + S. This algorithm is explicit and allows us to compute 
directly the points of each LS'-sequence. 

Some interesting extensions of LS-sequences of points to dimension two 
have been first introduced in [13] and presented in [I]. 

The sequence of points we want to study in this paper is associated to 
Kakutani's a-sequence {a n u} n >i, where a = ^p± is the inverse of the golden 
ratio. It is easy to see that, since a 2 + a = 1, each partition a n u is formed 
only by intervals of length a n (long) and a n+1 (short). 

If we denote by {l n } n >\ and {s n } n >i the number of long and short inter- 
vals, respectively, and by {t n }n>i the total number of intervals of the n-th 
partition, they satisfy the difference equation 

with initial conditions /o = 1 and l± = 1, so = and s± = 1 and, respectively, 
t — 1 and t\ = 2. Therefore {£ n }„>i is the Fibonacci sequence. For this 
reason, in [2] this sequence has been called the Kakutani-Fibonacci sequence 
of partitions, and it was also proved the following 

Theorem 5. The Kakutani-Fibonacci sequence of partitions {a n u;} n >i has 
low discrepancy. 

In [2] and [3] the first author proposed two algorithms which order in a 
natural way the left endpoints of {« n u;} n >i and associate to it a sequence 
of points called the Kakutani-Fibonacci sequence of points we will denote by 

{£n}n>l- 

The sequence {£ n }n>i can be constructed inductively as follows. 
Begin with the first block 

Ai = (6,60 = (0, a). 

If we denote by A n = . . . , £ tti ) the set of the t n points of a n u ordered in 
the appropriate way, the t n+ i — t n + l n points of a n+1 u are defined by 

An+l = (Ci--- ) Ct» ) ^n+l(Cl)>--- ) 0n+l(6„)) > W 
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I 1 1 1 1 1 1 1 1 

a 4 a 3 a 2 a 2 + a 4 a a + a 4 a + a 3 1 

Figure 1: First 8 points of the Kakutani-Fibonacci sequence {£ n }n>i 

where <p n+ i(x) = x + a n+1 . 

In Figure 1 we see the first 8 elements of the sequence {£ n }n>i- 
It is interesting to notice that the sequence {£„}„>i has been constructed 
by Ninomiya in [16J using a completely different approach, without observing 
that it is one of the possible reordering of the points of a Kakutani a-sequence. 

The next result has been proved independently by Ninomiya in [16J and 
by the first author in [2] using completely different methods. 

Theorem 6. The Kakutani-Fibonacci sequence of points {£ n }n>i has low 
discrepancy. 

In the next section we will obtain this sequence using a dynamical system 
approach. 

This approach to uniform distribution theory goes back to von Neumann 
and Kakutani (see [3], [TO], [IS])- F° r a recent survey on this matter see [TTj . 

Definition 7. Let X be a set and T : X — > X be a map. Given i6l, the 
sequence x, T(x), T 2 (x), ... is called the orbit of x in X. 

Definition 8. A map T : [0, 1[ — > [0, 1[ is said to be an interval exchange on 
[0, 1[ if there exist a finite or countable family J 7 of non empty subintervals 
h = [dk,bk[ of [0, 1[ and a family of corresponding real numbers Ck with 
1 < k < cardlJ 7 ), such that: 

(i) i h n h = and (4 + c h ) n (4 + c fc ) = if h ^ k, 

(ii) T(x) — x + Ck if x £ Ik, 

(iii) A(U fc >o4) = A(U fc > (4 + c k )) = I. 

It is easy to see that an interval exchange is a map of [0, 1[ into itself which 
is one-to-one, preserves the Lebesgue measure A and is continuous A-almost 
everywhere. 
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A special role is played by the ergodic transformations (see jH], for in- 
stance). 

Definition 9. Let (X,A,/-t) be a probability space. A measurable transfor- 
mation T : X — > X is said ergodic if for every A & A such that T~ l (A) = A, 
either /jl(A) = or = 1. 

We will use later an equivalent property: T is ergodic if and only if for 
any measurable set B such that /j,(B) > 0, the set B T = Ui=^oo T l {B) has 
measure 1. 

We will also assume that the ergodic transformation T is measure pre- 
serving, as most (but not all) authors do. 

The system (X, A, /i, T) is called a measure theoretical dynamical system, 
or dynamical system, for short. If T is ergodic, the system is called ergodic. 

The link between dynamical systems and uniform distribution is given by 
the following 

Theorem 10 (Birkhoff's Theorem). Let (X, A, /i, T) be a measure theoretical 
dynamical system. Then, for every f E C}{X) 

1 N ^ 
lim — y^f(T J x) 

j=0 

exists for fi-almost every x G X (here T°x = x). 
IfT-.X^-X is ergodic, then for every f G C l {X) we have 

l N ^ r 

lim - f(T j x) = / f{x)dn(x) , 

for fi- almost every x G X . 

The conclusion of the previous theorem says that the orbit {T^x}^ of 
x is u.d. for almost every x G X whenever T is ergodic. 

2 Main results 

We now turn to the main purpose of this paper showing that the Kakutani- 
Fibonacci sequence can be obtained as the orbit of the origin generated by 
an interval exchange T : [0,1[— > [0,1[, using the following cutting-stacking 
technique (see [9], [15] and [Tf]). 
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Definition 11 (Cutting- stacking technique for the Kakutani- Fibonacci se- 
quence). For every fixed n, the intervals of the n-th Kakutani-Fibonacci 
partition a n u are represented by two columns C n = {L n ,S n } constructed 
as follows. 

Let us start with two columns C\ = {L\,Si}, where L\ = [0,a[, Si = 
[a, 1[. If we divide L x proportionally to a and a 2 , we can write L\ = {L\, L\}, 
where L\ = [0,a 2 [ and L\ = [a 2 ,a[. Now we stack the interval Si over L\ 
(and use the common notation L\ * Si) as they have the same length, called 
width and denoted by w(L\) = w(Si) = a 2 . So we get two new columns 
C 2 = {L 2 , S 2 }, where L 2 = L°*Si and S 2 = L\. We denote by b(L 2 ) = [0, a 2 [ 
the bottom of L 2 and by h(L 2 ) — 2 — l 2 its height (or the number of intervals 
of which L 2 is made). 

L ° 2 

Ci 

L x Si " • " : ' S 2 



a 3 a 2 a 2 



v3 1 
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a a + a 3 a + a 3 1 



„ 4 



Figure 2: Partial graph of the cutting-stacking method for {a n u} n >i 
If we continue this way, at the n-th step we get two columns denoted by 



C n = {L n ,S n }, with h(L n ) = l n and h(S n ) = s n . We divide L n into two 
columns, say l? n and L\, where w(L° n ) = a n+1 and w{L^) = a n+2 , and stack 
S n over L°, obtaining 

Cn+l — {Ln+1, 'S'n+l}; 

where 
with 

w(L n+1 ) = a n+1 , h(L n+1 ) = l n+1 , b(L n+1 ) = [0, a n+1 [ 

and 

w(S' n+1 ) = a n+2 , h(S n+1 ) = s n+l , 6(C n+1 ) = [a n+1 , a n [ . 



The first steps of the procedure are visualized in Figure 2. 



To the above cutting-stacking construction it is naturally associated an 
interval exchange, whose explicit expression is given by the following result. 

Proposition 12. The interval exchange corresponding to the cutting- stacking 
procedure described in Definition [771 is the map T : [0, 1[— »■ [0, 1[ whose re- 
striction to I k is T k , where 



Ti(x) = x + a if x E h = [0, a 2 [ 



and, for every k > 1, 



k-l 



T 2k (x) = x + a 2k - a2j+1 if x G L 



2k 



k-l 



j=0 



j=0 



and 

T 2k +i{x) = x 



k-l 



+ a 2k+1 - « 2{J+1) if x £ I 2k+1 

3=0 



' k—1 k 

^a 2(j+1) ,^« 2 °' +1) 



3=0 



3=0 



Proof. If we write T(x) = x + c k whenever x e I k for any k > 1, we simply 
observe that I k + c k = [a k , a k + a k+1 [ for all k > 1. Therefore, A(Ufc>o(/fe + 
c k )) = 1 and (J/j + Cft) fl (I k + c k ) =0 whenever h ^ k, which proves that T 
is an interval exchange. 
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Now we show how the map T acts in the cutting-stacking procedure, 
proving that in each column C n = {L n , S n } the transformation T maps each 
interval of L n (respectively, S n ) onto the interval above it and the top interval 
of denoted as usual by top(L^), onto the bottom interval of S n . 

Let us start with C\ = {Li, Si}. When we divide the column L\ into two 
sub-columns L\ and L\, made by one interval each, we notice that top(L±) = 
Ii and b(Si) = T(I{). Therefore, we stack Si onto L® using T 1; which maps I\ 
onto b(Si) and gives birth to the columns C 2 = {L 2 , S2}, where L 2 = L®* Si 
and S 2 = L\. 

Now we consider n = 2k and C 2k = {L 2k , S 2k }. When we divide L 2k 
proportionally to a and a 2 , obtaining therefore L\ k and L\ k , we notice that 
b(S 2k ) = T 2k (I 2k ) = [a 2k , a 2k+1 [ and, consequently, top(L 2k ) = I 2k because 
T 2 k is a bijection. In other words, T 2k stacks the bottom of 5*2^ onto the top 
of L° 2k because T 2k (top(L° 2k )) = b(S 2k ) = [a 2k , a 2k ~ l [ . 

A simple calculation shows that if we consider the case n = 2k + 1 we have 
T 2k+ i(I 2k+ i) = b(S 2k+ i) = [a 2k+1 ,a 2k [ and, therefore, top(L° 2k+1 ) = I 2k +i. 

The proposition is now completely proved. □ 

Definition 13. The map T : [0, 1[— > [0, 1[, whose expression is given in 
Proposition [121 is called the Kakutani-Fibonacci transformation. 

Figure 3 shows the graph of the maps T k , with 1 < k < 4. 

Note that in [TTJ it is presented a transformation Tf, called the Fibonacci 
transformation, related to the translation x — > x + a(mod 1) by means of the 
cutting-stacking technique. At each step of the cutting-stacking procedure, 
the two columns obtained by means of Tf have the same height and width 
as the two columns obtained by means of the Kakutani-Fibonacci transfor- 
mation T, but the bottoms are different as well as the order of the intervals 
in each column. Moreover, Tf is not defined in 0. 

Theorem 14. The Kakutani-Fibonacci transformation T is ergodic. 

Proof. The arguments we use are inspired by [9]. 

Let us denote by C* m the stack made by L m * S m and let us define the 
mapping r m : [0, 1[— > [0, 1[ as follows. For 1 < i < l m — 1 it is the translation 
of the interval Ji of L m onto Ji+i, the interval just above Jj, and hence it 
coincides with T. If i = l m , it is a contraction by parameter a of top(L m ) 
onto b(S m ). For 1 < i < s m — 1 it is the translation of the interval Jj of S m 
onto Jj+i, the interval just above Jj, and hence it coincides with T again. 
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Figure 3: Partial graph of the Kakutani-Fibonacci transformation T 

Now we fix a measurable set B in [0, 1[ with X(B) > 0, such that T^ 1 {B) = 
B, and we prove that X(B T ) = 1, where B T = Ut-oo T\B). 

From Lebesgue density theorem, there exists a point xq having density 1 
for B. This implies that for every fixed e > there exists 5 > such that if 
I is any interval containing x with X(I) < 5, we have 

X(B n J) > (1 — e)A(J). (2) 
Due to the fact that the diameter of the Kakutani-Fibonacci sequence of 
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partition (a m u) m >i tends to when m — > oo, there exists m > such that 
a m < 5 and, therefore, there exists an interval of for which ([2]) holds. For 
sake of brevity, we will denote it by /. 

As = L m * S m , the interval I could belong to the column L m or to 
the column S m . If / £ Sm, it is equivalent to say that I G L m+ \. For this 
reason, without any loss of generality we can suppose that / £ L m . 

Suppose now that / is the r-th interval from below in the column L m and 
observe that 



t n 



M U r ™( J ) = E A ( J ) = 1 - 

\i=— r+1 / i=— r+1 

Taking the above identity and ([2]) into account, we have the following 
inequalities: 

X(B T ) > \{(Bni) T ) 

(hm—r \ Am-1 

i=l-r / V i=0 

(t m -r \ / im-r 

U 4(£n/) +A |J r^(sn/) 
i=l-r / \7={m+l-r 

= E A ( r -( Sn/ ))+ E A K( 5 n/)) 

i=l— r j=l m +l—r 

> (l-e)l m a m + (l-e)s m a m+1 = l-e. 

As e is arbitrary, we conclude that X(B) = 1. Therefore, T is ergodic and 
the theorem is proved. □ 

A direct consequence of the above theorem and Theorem [10] is the follow- 
ing 

Theorem 15. The sequence {T ra (x)} n >o is u.d. for almost every x G [0, 1[. 

We conclude this paper showing how to get the points of the Kakutani- 
Fibonacci sequence by means of the Kakutani-Fibonacci transformation. 

Theorem 16. The Kakutani-Fibonacci sequence of points {^ n }n>i coincides 
with {T"(0)} n > . 
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Proof. We want to show by induction on n > 1 that 

(6, 6,..., et B ) = (o,r(o),T 2 (o),...,r t »- 1 (o)). (3) 

If = 1, Q is obviously verified. 

We suppose that (j3J) is true for n and prove that 

( £l, $2, 6„> 0n+l(£l)> • • ■ , 0n+l(6n)) 

(0, T(0), T 2 (0), . . . , T*" _1 (0), T tn (0), T* n+1 (0), . . . , T^+^^O)) . 
Due to ([I]) and the inductive assumption, it is sufficient to prove that 
(0„+i(O), 0n+i(T(O)), . . . , ^(T^-^O))) 

(r* re (o), r*' i+1 (o), . . . , r tn+,n - 1 (o)) . 

As <j) n+ i(x) = x + a n+1 , the above identity is equivalent to 
(a {n+1) , T(0) + a {n+l \ . . . , T^-^O) + « (n+1) ) 

= (T* n (0),T* n+1 (0),...,T* n+l -- 1 (0)). (4) 

In order to prove 01]), we focus our attention on the intervals of the column 
S n+ i and on their left endpoints. 

We note that, due to the cutting-stacking procedure described in Defi- 
nition [TT] and to the nature of T described in the proof of Proposition [12] and 
specifically to the fact that b(S n+1 ) = T(top(L n+1 )) = T(T / " +s ™- 1 ( [0, a n+l [ )), 
the columns C n+ i = {L n+ i, S n+ i} can be written as follows: 

L n+1 =L° n *S n = (b(L° n ) * ■ ■ ■ * top(L°Sj * (b(S n ) * ■ ■ • * top(S n ) 
= ([0,a n+1 [ *•••* T ln ~ 1 ( [0,a n+1 [ ) ) * 
* (T ln ( \0,a n+1 [) *-•• * T l - +s --\ [0,a n+1 [)) 

and 

S n+1 — L\ = b{L l n ) * ■ • • * tap{Ll) 

= (\a n+ \ a n [ * ■ ■ ■ * T 1 "- 1 ( [a n+ \ a n [ ) ) 
= [0,a n+1 [ ) * ••• * T l n+in-i ( [ 0ja «+i[ ) y 
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Therefore, as T is an interval exchange, the t n left endpoints of L n+1 
are 0,T(0),--- , T* n_1 (0) and the left endpoints of S n+ i, which are a right 
shift by the constant a n+1 of the left endpoints of L n+1 , are T* n (0), T* n+1 (0), 
■ ■ • ,T* n+s "~ 1 (0), which proves fl3J). 

The theorem is now completely proved. □ 

Acknowledgements The authors wish to express their thanks to Robert 
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